Order parameter node removal in the d-wave superconductor Y Ba2Cus07-x under 

magnetic field 
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Whether the node in the order parameter characteristic of a d — wave superconductor can or 
cannot be removed by an applied magnetic field has been a subject of debate in recent years. 
Thermal conductivity results on the high Tc superconductor Bi2Sr2C aCu20s originally explained 
by Laughlin in terms of such a node removal were complicated by hysteresis effects, and judged 
inconclusive. We present new tunneling data on Y Ba2Cus,07-x that support the existence of the 
node removal effect, under specific orientations of the sample's surfaces and magnetic field. We 
also explain the hysteretic behavior and other previous tunneling results so far not understood 
satisfactorily, attributing them to a combination of node removal and Doppler shift of low energy 
surface bound states. 

PACS numbers: 74.50.+r, 74.72.Bk 



The excitation spectrum of a conventional supercon- 
ductor (Low Tc) is characterized by an almost indepen- 
dent momentum (s — wave) energy gap, A. This is not 
the case in the high Tc cuprate superconductors; there 
it is broadly agreed that the ground state superconduct- 
ing order parameter (OP) is strongly momentum depen- 
dent, it is maximal in the direction of the crystallographic 
axes a and b; in most cases it appears to have the pure 
dx2-y2 symmetry, being zero at 45° between these axes 
(the node directions), where it changes sign. In contrast 
to the finite energy A required to excite a low energy 
quasi-particle in a Low Tc superconductor, such a quasi- 
particle can be excited with an infinitely small energy in 
ad — wave superconductor along the nodes. This is no 
longer the case if an additional imaginary component is 
present in the OP, on that case the energy spectrum of 
the superconductor is fully gapped. 

Theoretically, it has been suggested that such an 
imaginary component can result from instability of the 
d — wave OP under perturbations such as surface pair 
breaking impurities |2j, proximity effect Q and 
magnetic field Another view is that a phase transi- 
tion occurs at a certain doping level or magnetic 
field Q from pure d — wave to a nodeless OP having the 
dx2-y2 -\- idxy or d.j;2_y2 -\- is symmetry. The idxy com- 
ponent breaks both time and parity symmetries, hence, 
as pointed out by Laughlin |^ , it involves boundary cur- 
rents that flow in opposite directions on opposite faces 
of the sample. They produce a magnetic moment which, 
through interaction with the magnetic field, lowers the 
free energy by a term proportional to B ■ dxy, where B is 
the magnetic field. On the other hand, in the zero tem- 
perature limit, node removal costs an energy proportional 
to Minimization of the sum of the two contribu- 
tions leads to an amplitude \dxy 
a coefficient. 



= A ■ B^/^, where A is 



Experimentally, two sets of experiments have been in- 



terpreted as indicating that a magnetic field, applied per- 
pendicular to the Cu02 planes, can indeed induce a node 
less OP. Measurements of the thermal conductivity k{H) 
on Bi2Sr2CaCu20s (Bi2212) single crystals have shown 
a decrease followed by a plateau at a certain field 
This field was interpreted as that beyond which a finite 
idxy component appears at the finite temperature where 
the experiment is performed ||^. The finite gap in the 
plateau region prevents the excitation of additional quasi 
particles. Field hysteresis of the thermal conductivity has 
however led to a controversy as to the actual origin of the 
plateau, and this issue has remained unresolved ||lO|| . 

A second set of experiments possibly indicating the 
occurrence of a nodeless OP is the field evolution of the 
conductance of in-plane tunnel junctions formed at the 
surface of YBa2Cu307^x (YBCO) films oriented per- 
pendicular to the Cu02 planes |23| . This conductance 
presents a peak at zero bias (known as Zero Bias Con- 
ductance Peak, ZBCP), which splits when the field is 
applied parallel to the film's surface, and perpendicular 
to the Cu02 planes. The ZBCP is one of the clearest 
manifestations of an OP having nodes, because it comes 
about due to a change of phase by tt upon reflection at a 
(110) surface, which generates low energy surface states, 
or Andreev bound states (the fact that ZBCPs are 
also observed for other in-plane orientations is generally 
interpreted as due to surface roughness ||ll| ) . Peak split- 
ting may indicate node removal, the new peak position, 
S{H), giving the amplitude of the dxy component 
However, a different explanation of the peak splitting 
was also offered, in terms of a Doppler shift of the An- 
dreev bound states This Doppler shift is equal to 
Vg-pp., where Vg is the superfluid velocity associated with 
the Meissner currents, and the Fermi momentum of 
the tunneling quasi-particles. There are difficulties with 
both interpretations. The amplitude of the idxy compo- 
nent should be essentially reversible at the fields of inter- 
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est (in the Tesla range), where the thin film samples are 
well into the Bean complete penetration limit, while in 
fact a strong field hysteresis is observed |12, llM H. As 



for the Doppler shift, it should vanish at film thickness 
smaller than the London penetration depth where the 
Meissner superfluid velocity is much reduced, varying as 
Vs = eXHtanh{d/2X), where d is the film thickness, A 
the penetration depth and H is the magnetic field at 
the sample surface; while in fact it is slightly changes at 
d = (A/2) |2^. A further difficulty is that, for reasons 
that have remained unclear until now, field splitting of 
the ZBCP is not always observed. For instance, it was 
not seen in Lai.85<S'ro.i5CM04 and YBCO grain bound- 
ary junctions [15l, nor in junctions prepared on Bi2212 
single crystals [ [f6[ . Thus, a consensus has not yet been 
reached as to whether a field can induce a finite sub-gap 
in a d-wave superconductor. 

We present in this Letter a new series of tunneling ex- 
periments that clarify the origins of the field splitting 
phenomenon, and establish the conditions under which 
node removal occurs. Our central result is that node re- 
moval can be observed in (110) oriented samples only. 
It is most clearly seen in decreasing fields, for which the 
ZBCP splitting is not affected by Doppler shift effects, 
and is thickness independent. Doppler shift does affect 
data taken in increasing fields, which are thickness de- 
pendent. Hysteresis is due to the properties of the Bean 
Livingston barrier, which is effective against flux pene- 
tration in increasing fields, but not against fiux exit in 
decreasing fields. Data taken on (110) oriented samples 
in decreasing fields, starting from a field of up to 16T, 
are in quantitative agreement with Laughlin theory. 

YBCO films nearly optimally doped, having thickness 
ranging from 600Aup to 3200A, were prepared in the 
(110) and (100) orientations by the template method, 
using SrTiOs and LaSrGaO^ substrates of the appropri- 
ate orientation ||l^, |l8j . Critical temperatures of all films 
were in the range of 88K to 90K. Junctions were prepared 
by pressing In (Indium) pads on the films fresh surface 
0, |18|. AU junctions were measured at 4.2K, and some 
also at 1.6K. All junctions displayed an unsplit ZBCP in 
zero magnetic field, irrespective of the film orientation. 
Junctions characterisitics were measured as a function 
of field, applied parallel to the surface, either parallel or 
perpendicular to the Cu02 planes. A typical data set 
is shown Fig.l for a 600A, (110) oriented film. Mea- 
surements were systematically taken in increasing and 
decreasing fields. Field splitting was only observed when 
the field was applied perpendicular to the Cu02 planes, 
in agreement with previous results, confirming the uni- 
axial in-plane orientation of the c-axis |l^ . A total of 20 
junctions were measured. 

The peak position S{H) measured on the 600 A thick, 
(110) oriented film whose conductance characteristics are 
shown Fig.l, is plotted Fig. 2a in increasing and decreas- 
ing fields. In increasing fields, the peak position reaches 



about 4 meV at 5T. At higher fields, it cannot be de- 
termined because the peak becomes too smeared, possi- 
bly because it merges with the main gap structure. For 
thicker films, the rise in the peak position is faster. For 
a 3000A thick film, it reaches 4 meV already at 3T. Data 
taken in decreasing fields for the 600A film starts at about 
4 meV at 15T. By contrast with the behavior in increas- 
ing fields, that in decreasing ones is thickness indepen- 
dent. 

The hysteresis amplitude is shown Fig. 2b up to 5T, the 
highest field at which it could be determined. Hysteresis 
saturates above 2T, a field of the order of He. For thicker 
films, the hysteresis amplitude reaches larger values, but 
the field dependence is similar. This hysteretic behav- 
ior can be understood within the Doppler shift model 
of the ZBCP splitting if we take into account the 
properties of the Bean Livingston barrier. This barrier 
prevents flux entry up to fields of the order of He. Up 
to that field, the superfiuid velocity of the Meissner cur- 
rents increases almost linearly; beyond that field, it sat- 
urates. The Doppler shift, proportional to the superfiuid 
velocity, follows the same behavior. By contrast, in de- 
creasing fields, there is no barrier that prevents flux exit. 
This has been shown theoretically by Clem , verified 
experimentally by Bussiere | pO{ , and shown also to ap- 
ply in tunneling experiments by Moore and Beasley [ pl| . 
Hence, when the field is decreased, the surface superfluid 
velocity quickly reduces to zero, and so does the Doppler 
shift. This is the origin of the hysteresis. The fact that 
the hysteresis amplitude increases with fllm thickness is 
in line with this interpretation. 

Our central result is shown Fig. 3, which presents the 
ZBCP splitting measured in decreasing fields for (110) 
oriented films having thickness ranging from 3200A down 
to 600A, plotted as a function of the square root of 
the applied field. Data for all samples follow the law 
Si ^ A- B^/^ with A = l.lmey/ri/2. Laughlin calcu- 
lated for Bi 2212 a coefficient A = 1.6meV/T^/'^. This 
coefficient is proportional to the square root of the gap 
and that of the Fermi velocity. Taking into account that 
gap values measured on Bi2212 (25 to 30 meV) are larger 
than those measured on YBCO (17 to 20 meV), and as- 
suming that values of the Fermi velocity are similar for 
both compounds near optimum doping, there is a good 
quantitative agreement between theory and experiment. 

These results are quite different from those previ- 
ously reported for (100) |l§ and (103) |l|, H, ^ ori- 
ented films. Generally speaking, field splitting values are 
smaller for these orientations. For 3000A (100) oriented 
films, the splitting measured in increasing fields reaches 
only 1.5 meV at 5T Q, as compared to 4 meV at 3T for 
a similar thickness in the (110) orientation, as mentioned 
above. Differences are even larger in decreasing fields, 
splitting values falling always well below the line shown 
in Fig. 3. For the 3000A thick sample mentioned above, 
it is only of 0.5 meV at 5T, as compared to 2.4 meV in 
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the (110) orientation. In addition, splitting values in the 
(100) orientation are strongly thickness dependent. For 
samples thinner than I6OOA, we find that they are too 
small to be determined experimentally at 4.2K up to 6T. 
These results are consistent with a splitting dominated 
in the (100) orientation by the Doppler shift effect, with 
the Laughlin mechanism playing apparently no role. 

Surface faceting is thought to be the reason for the 
ZBCP commonly observed in (100) oriented films 
Recently, some direct evidence has indeed been provided 
by STM measurements, suggesting that (110) faces are 
present in films of that orientation ||l^. While the zero 
field ZBCP is similar for both macroscopic orientations, 
their field splitting is entirely different, as reported here. 
Our results demonstrate that in order to observe a sub- 
stantial splitting in decreasing fields and in films smaller 
than the London penetration depth - namely, under con- 
ditions where the Doppler shift effect is very weak - one 
must use samples having the (110) orientation. Then, 
and only then, the experimental data are in agreement 
with the law of Laughlin, strongly suggesting that node 
removal does occur. Even though (100) oriented films do 
have (110) facets at their outer surface, the interface with 
the (100) LaSrGaO substrate and PrBaCuO intermedi- 
ate layer is presumably quite fiat. We conjecture that 
the absence of the second (110) surface prevents the flow 
of Laughlin's currents on that face, making the establish- 
ment of the idxy component energetically unfavorable. 

In conclusion, the presented tunneling data is consis- 
tent with node removal in the d — wave superconductor 
YBCO under magnetic field if, and only if, the sample's 
boundaries have the (110) orientation. As far as we know, 
this has not been predicted theoretically. Previously not 
well understood experimental results, such as the differ- 
ence in field splitting between (110) and (100) oriented 
films and its hysteretic behavior, can now be explained. 
The absence of the ZBCP field splitting in grain bound- 
ary junctions IT^ can be understood as a combination of 
two factors: first, the field being applied perpendicular 
to the surface, vortices penetrate at low fields, and there 
cannot be any substantial Doppler shift. Second, the ge- 
ometry of the boundaries is unfavorable for the flow of 
Laughlin's currents. The same applies to junctions pro- 
duced at the edge of Bi2212 crystals jl^ with the field 
being applied perpendicular to the surface. It could be 
that the contradictory results reported in thermal con- 
ductivity experiments, concerning the existence of a field 
induced gap, also stem from the ability or inability of 
the samples to carry boundary currents under the spe- 
cific experimental conditions. 
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FIG. 1: Normalized dynamical conductance G=dI/dV Vs bias 
V for increasing (a) and decreasing (b) applied magnetic fields 
for YBCO (110) oriented film at 4.2K. Film characteristics: 
Tc= 88K, film thickness d=600A. The splitting (5) is defined 
as half of the distance between the positions of the conduc- 
tance maxima. In increasing field it can be determined clearly 
from field of about O.IT up to 5T, and in decreasing fields 
from 13T down to 0.9T. Applied fields (in Tesla) : 0, 0.1, 0.3, 
0.5, 0.7, 0.9, 1.2, 1.5, 1.8, 2.1, 2.5, 3.0, 3.5, 5, 6, 7, 11, 13, 15. 
(c) behavior of the same junction for magnetic field applied 
parallel to the CuO planes at fields (in Tesla): 0, 0.5, 1, 2, 4, 
8, 12 ,15.5. The strong anisotropy of the field effect confirms 
the good in-plane orientation of the c-axis. 
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FIG. 2: Comparison between the field splitting hysteresis 
curves for (110) and (100) in the presence of increasing (A) 
and decreasing (v) external magnetic fields, (a) ZBCP split- 
ting (5) and (b) ZBCP splitting difference {A6 = 6^ - di) for 
(110) 600A thickness films at 4.2K. The line in (b) is a guide 
to the eye. 
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FIG. 3: ZBCP Field splitting measured in decreasing fields 
for film thickness ranging from 3,200A to 600A as a function 
faceting square root of applied magnetic field (H). The line is 
a linear fit to all points with l.lmV/T^^^ slope. 



